1/Ti in the D-wave superconducting state with coexisting antiferromagnetism 
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We consider a general D-wave superconducting gap function in the presence of coexisting antifer- 
romagnetic (AFM) long-rang order. We find that the D-wave order parameter develops additional 
nodes without lowering the symmetry of the state due to the interplay of the superconducting and 
AFM correlations. This AFM correlated D-wave gap has a small gap structure inside the generic 
D-wave gap. As a result of this additional structure it shows a quite different response to impurities 
compared to a standard D-wave gap. We calculate the nuclear spin-lattice relaxation rate 1/Ti of 
this gap in the presence of impurities and discuss the implications for the experiments of some AFM 
heavy- fermion superconductors and in particular for CeRhIn5. 
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I. INTRODUCTION 



The interplay and competition between antiferromag- 
netism (AFM) and superconductivity (SC) is a long 
studied subject both experimentally 1 and theoretically 2 . 
There are now several heavy fermion compounds show- 
ing clear experimental evidence of the microscopically ho- 
mogeneous coexistence of AFM and SC correlations 3-6 . 
However, the direct consequences of the coexistence of 
AFM and SC has not yet been clearly resolved. Re- 
cent 115 In- NQR studies of CeRhlns and Celns also re- 
vealed that there is a coexistence region of AFM and SC 
around the critical pressures of 1.6 GPa and 2.43 GPa, 
respectively 7 . Apart from confirming the coexistence of 
AFM and SC in this region of pressure, an interesting ob- 
servation is that the nuclear spin-lattice relaxation rate 
1 jT\ displays the linear temperature dependence 1 jT\ oc 
T at low temperatures (up to 0.3 T sc for CeRhlns), fol- 
lowed by the T 3 behavior below the superconducting 
transition temperature T sc . The latter behavior is an in- 
dication of the lines of nodes of the anisotropic SC state. 

The standard explanation of the T-lincar temperature 
dependence of 1/Ti at low temperatures is the impu- 
rity effect in an unconventional SC with lines of nodes. 
However, this impurity explanation has difficulties in the 
case of CeRhIn 5 ; the low temperature T- linear 1/Ti is 
only observed in the AFM and SC coexistence region at 
1.6 GPa, and once the AFM order disappears by increas- 
ing pressure to 2.1 GPa, 1/Ti exhibits a T 3 behavior for 
the entire measured temperatures below T sc 8 . Similar 
pressure dependence of 1/Ti has also been observed in 
Ge-substituted CeCu2Si2 5 . Another AFM heavy fermion 
superconductor URu2Si2 4 also shows the low tempera- 
ture T-linear behavior in 1/Ti. In the case of UPd2Al3, 
the situation is somewhat different; a more recent NMR 
experiment 6 indicates that the T-lincar behavior ob- 
served in earlier experiments disappears with purer sam- 
ples. Putting together these observations we think that 
the interplay of the AFM and SC, together with some 
sort of impurity effects, should hold the key to explain 



the 1/Ti oc T behavior at low temperatures in these AFM 
heavy fermion superconductors. In particular, we think 
this is the case for CeRhIn 5 under pressure, where at low 
temperatures 1/Ti oc T is observed in the coexistence 
region (P=1.6 GPa), but is replaced by 1/Ti oc T 3 once 
the AFM is suppressed (P=2.1GPa). 

There is a substantial body of studies on the inter- 
play/coexistence between D-wave SC and AFM 2 . The 
general conclusion of these studies is that, depending on 
the parameters in the model, these two phases either ex- 
clude each other or can coexist together. Here, how- 
ever, we are motivated by the above-mentioned magnetic 
heavy-fermion (HF) superconducting materials, where 
the Neel temperature TV usually sets in first at a higher 
temperature and then at a lower temperature the SC 
transition occurs on top of the AFM ordering. Hence- 
forth, we assume the coexistence of D-wave SC and AFM 
and confine our interest on how a general D-wave SC or- 
der parameter (OP) should be modified due to the pres- 
ence of the AFM long range order. We found an impor- 
tant modification of the D-wave gap due to the long range 
AFM order; the D-wave SC state develops additional 
nodes in the gap function that are a direct consequence 
of the crossing between the Fermi surface (FS) and the 
magnetic Brillouin zone (BZ), see Fig.l. At these dis- 
crete crossing points the superconducting order param- 
eter is forced to vanish by the AFM coherence factors 9 . 
These additional nodes (we call them AFM nodal points) 
lead to a density of states (DOS), which is almost gap- 
less down to very low energy before it develops a steeply 
decreasing DOS, see Fig. 2. This shape of the DOS is 
intrinsically vulnerable to a small amount of impurities, 
which fill rapidly the small gap and produce a large finite 
DOS. 

From this AFM+D-wave gap solution we calculate the 
nuclear spin-lattice relaxation rate 1/Ti in the presence 
of impurities. We find that the response of this AFM+D- 
wave gap to impurities is quite different from the one of 
the standard D-wave gap. The main findings are: (1) the 
AFM+D-wave gap produces more quickly a finite DOS at 
oj = than the standard D-wave gap, both for Born and 
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for unitary impurity scattering, (2) for unitary scatter- 
ing, the calculated 1/Ti, both for the AFM+D-wave gap 
and for the standard D-wave gap, produces the T-linear 
behavior at very low temperatures. The temperature re- 
gions of T-linear behavior are about the same for both 
gaps with the same amount of impurity concentration 
despite the different low energy behaviors of the DOS. 
However, the magnitude of 1/Ti in the T-linear region is 
about 4 times smaller for the standard D-wave gap than 
for the AFM+D-wave gap, (3) for Born limit scattering, 
the l/Ti of the AFM+D-wave gap can still produce the 
T-linear behavior in a substantial region of low temper- 
atures but with a much higher impurity concentration 
than with unitary impurities. In order to achieve the 
same T-linear behavior as for unitary scattering, about 5 
times the amount of impurity concentration is needed. 
On the other hand, the standard D-wave gap hardly 
shows any T-linear behavior with the same amount of 
impurities, but displays almost T 3 power law for the en- 
tire temperature range below T sc with a small deviation 
due to impurities. 

Considering the experimental observation of the T- 
linear behavior in 1/Ti in some of the AFM heavy 
fermion superconductors, we propose that the AFM+D- 
wave gap structure together with Born limit impurities 
can be the reason for the T-linear behavior in 1/Ti in 
these compounds. When the AFM long range order dis- 
appears by tuning pressure in CeRhlns, as well as in 
CeCu2(Sio.98Geo.o2)2, the SC gap gradually returns to a 
standard D-wave gap and then a small amount of Born 
scatterers is not enough to produce the T-linear 1/Ti 5 ' 7 . 
One remaining question of this scenario is the possible 
source of the Born limit scatterers. We think that the 
most probable source of scatterers in the Born limit is 
strain in the lattice induced with pressure and/or AFM 
domain walls in CeRhlns and through chemical substi- 
tution in CeCu 2 (Sio.98Ge .o2)2 5 - 



II. FORMALISM 

The Hamiltonian describing the coexistence of AFM 
and SC pairing interactions is written as 

H = ^2,[e k c\ k(y c k ^ + aA M c} k+Q<J c k ^} 

+ V{ k i k')tf k ,aC\-k,a>C-k>,a>C k >^, (1) 
k,k' 

where we replaced the AFM pairing interaction by its 
mean- field AFM OP A M . In the presence of the AFM 
order, the conduction band is split into two AFM quasi- 
particle bands: 

a (l),a,k = u k c a _ k — av k c a ^ k+ Q, 

a (2),<7,k = WkC^k + ItfcC^fc+Q, (2) 



with the AFM Bogoliubov coefficients, u\ = 1/2 + 

, j h , , and vl = 1/2 , t* , ■ The quasiparti- 

cle energy and the convenient parameters are defined as 

E(1 U 2) = Cfc ± V^ + A m> (3) 
Ck = ^[efc + e k+Q ], (4) 
Sk = - e k+Q }. (5) 

Then the original SC paring interaction is rewritten in 
terms of the AFM quasiparticle wave function, 

H in t = ^2v{k,k')c-\ 1k c-\ l _ k c l - k >c hk > 

k,k' 

= V ( k > k ')( U l - W fc) a t(2)T,fc Q t(2)l,-fc 

k,k' 

a(2)i,-fc'a(2)T,fc'( u fe' - v l>) + •••> ( 6 ) 

where the ellipsis stands for terms like o<xt Q! iti Q! i Q; i an d 
all the mixed terms of a\ and ai. In this approximation 
we certainly keep the pairing term in the lower AFM 
quasiparticle band assuming that the Fermi level crosses 
the lower band and that the AFM gap Am is much larger 
than the SC gap A sc . This approximation should apply 
for the case TV 3> T sc . Then the singlet gap equation is 
written as 

A a (k) = (u%\ — v%) V(k, k') < a (2)hk ,a { 2)i,-k> > 
w 

Aul-vl). (7) 

This gap equation was obtained by Buzdin and 
Bulaevskii 9 for s-wave SC. The important point is that 
the AFM coherence factor (u\ — vf.) modifies the other- 
wise standard D-wave gap equation where SC is medi- 
ated by the pairing potential V(k,k'). At those points 
where {u\ — vj:) = is satisfied, the gap function A a (k) 
should vanish in addition to the D-wave nodal points. 
Now for simplicity of numerical calculations, we assume 
a circular FS in two dimensions and the gap equation 
is generalized to account for impurity effects 10 . The ef- 
fect of the impurity scattering is included with T-matrix 
approximation 11 . For the particle-hole symmetric case 
Tj = 0, and for D-wave OP with isotropic scattering T\ = 
(also without loss of generality we can choose T 2 = by 
U(l) symmetry). Then we need to calculate only Tq(u>). 
The impurity selfenergy is given by Eo = TTo, where 
r = rii/irNo, Nq is the normal DOS at the Fermi energy, 
rii is the impurity concentration; To(uj n ) = ^ i 

where g Q {u> n ) = ^ E fc j+A^(fc) ■ Q n - + S 
(u> n = 7r T (2n+ 1)), and the scattering strength param- 
eter c is related to the s-wave phase shift 6 by c = cot (J). 
With this To the following gap equation is solved self- 
consistently. 
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A(0) = -N(0) ■ g{4>) I -^V{<p-<i>) 



) 



l(4>), (8) 



numerical calculation (3 = 1.74 is fixed because this pa- 
rameter is not sensitive for the final results. But the ratio 
Ao /T sc is an important parameter to simulate strong cou- 
pling effects; the larger this ratio is the more the strong 
coupling effect is accounted for. 



where g(<j>) is the angular parametrization of (u 2 . — v 2 ,). 
In contrast to previous works on the D-wave SC, where 
the form of the D-wave SC gap has a fixed functional 
form such as A(k) = A (cosfc 2: — cosk y ) or cos(20), we 
allow for a most general D-wave OP of T>2 symmetry; 
namely A(mr/4) = 0(n = 1,3,5,7), A(</>) = A(0±tt), 
and A(0) = — A(<fi ± tt/2). Therefore the gap equation 
can produce the most general D-wave symmetry gap so- 
lution for a given pairing potential and now with an addi- 
tional constraint from the AFM correlation. The pairing 
potential V(<f> — <f> ) induces a D-wave gap. Although 
its microscopic origin is not an issue in this paper, we 
believe it originates from AFM fluctuations. The static 
limit of AFM fluctuations x(q, u> = 0) ~ ^ q _Qyi + ^--i is 
parameterized as 10,12 



V(<t>-<f>) = V d (b)- 



±n/2) 2 + b 2 



(9) 



where the parameter b is proportional to normalized 
in the circular FS (£ ~ air/b; a is the lattice parameter). 
For all calculations in this paper, we choose b = 0.5 which 
is not a sensitive parameter for our results. 

With the gap function A(<j>) and To(u>) obtained from 
Eq.(8) (7o(cj) is analytically continued from T (uj n ) by 
Pade approximant method), we calculate the 1/Ti nu- 
clear spin-lattice relaxation rate 11 ' 13 
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(10) 



where Co = co + Y,$(lo) and (...)^ means the angular av- 
erage over the FS. The first term in the bracket of Eq. 
(10) is N 2 (io). The second term vanishes in our calcu- 
lations because of the symmetry of the OP. To calcu- 
late \jT\T using Eq. (10), we need the full temper- 
ature dependent gap function A(<fi,T) and T sc . Our 
gap equation Eq. (8) is basically the BCS gap equa- 
tion, therefore it gives the BCS temperature behavior 
for A((f>,T) and A /T sc = 2.14 for the standard D-wave 
SC. In order to account for strong coupling effects and 
a stronger anisotropy of our general D-wave gap solu- 
tion, we use the phcnomcnological formula, A(</>, T) = 
A(<p, T = 0) E(T) with S(T) = tanh(/3^T sc /T - 1), and 
parameters j3 and A /T sc . Then we only need to calcu- 
late A((j>, 0) at zero temperature. The temperature de- 
pendence of S (w, T) (= r Tq(oj, T)) is similarly extrapo- 
lated: T (u,T) = T (lo,T = 0) S(T)+r normoi (l-S(T)), 
where T norma i = T/(c 2 + 1) is the normal state T . In our 



III. RESULTS AND DISCUSSIONS 

Now let us discuss the numerical results. In Fig. (1) 
the schematic FS and the magnetic BZ are shown to- 
gether with the two-dimensional original BZ. The points 
where the magnetic BZ and the original FS cross and 
therefore the condition (u 2 — v 2 ) = is satisfied, are 
marked as AFM nodal points. The gap function vanishes 
at those points; in real space this means that the singlet 
SC pairing is prohibited along the (1,1) direction due to 
AFM spin correlations 9 . A schematic D-wave gap solu- 
tion is drawn accordingly. 

In Fig. (2), we show the normalized DOS N(uj)/Nq for 
two exemplary cases of the AFM nodal points (77 = 0.2 
and 0.4); the distance of the AFM nodal point from the 
D-wave nodal point is parameterized by 77 = \4>afm ~ 
4>d\/4>d, where 4>d = 7r/4 is the D-wave nodal point and 
4>afm is the AFM nodal point. The inset shows the 
corresponding gap solutions. Wc found empirically that 
A((/>) ~ [cos(20) + 1.5 77COs(6</))] is a very good approx- 
imation to fit the AFM+D-wave gap solutions. In real 
compounds, the distance of the AFM nodal points from 
the D-wave nodal points is determined by the shape of 
the original FS (before the AFM long range order is in- 
cluded). The key feature is that when the AFM nodal 
points are not very far from the D-wave nodal points, 
the DOS already develops the gapless feature (without 
impurities), except at very low energies. This very low 
energy region has a shallow gap (the steeply rising DOS 
region), which is intrinsically vulnerable to disorder. 

Fig. 3 (a) shows the DOS and 1/Ti calculated for the 
AFM+D-wave gap with unitary impurities (c = 0) . The 
inset shows the normalized DOS of the AFM+D-wave 
gap solution with varying concentrations of the unitary 
scatterers. We choose 77 = 0.3 (see Fig. 2) for illustration 
of a typical AFM+D-wave gap. As explained above, the 
low energy region of the DOS is quickly filled with a small 
amount of impurities. With this choice of parameters 
the impurity scattering rate r/A = 0.064 is enough to 
completely fill the low energy gap region and N(lo = 0) 
reaches more than 50 % of the normal-state DOS No. 
The main panel shows the nuclear spin-lattice relaxation 
rate 1/Ti in a log-log scale for the corresponding DOS in 
the inset. In Fig. 3 and Fig. 4 all 1/Ti results are nor- 
malized to 1/Ti = 10 at T=T SC for easy comparison. For 
the temperature dependence of the gap, we choose the pa- 
rameters j3 = 1.74 and A /T sc = 3 for the AFM+D-wave 
gap to account for the strong coupling effects of super- 
conductivity as explained before. As expected from the 
DOS results, due to the impurity induced residual states, 
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l/Ti displays the lincar-T dependence at low tempera- 
tures and the region of T-linear behavior increases with 
impurity concentration. For T/Aq = 0.064. this T-lincar 
region extends up to ~ 0.35 T sc . In the high temperature 
region, first of all, near T sc the coherence peak is almost 
invisible because of the sign-changing gap function (the 
second term in Eq. (10) vanishes). Secondly, it shows the 
typical T 3 behavior below T sc due to the lines of nodes in 
the gap until it goes through a gradual crossover region 
and finally to the T-linear region. The comparison with 
the experimental data of CeRhlns at 1.6 GPa 7 (green di- 
amond symbols; the experimental data are normalized in 
the same fashion as the numerical results) is quite good 
with the theoretical calculation with T/Aq = 0.064. On 
the other hand, there is no way to fit the 2.1 GPa data 
(magenta diamonds) with any amount of impurities. 

Fig. 3 (b) shows the same calculations as in Fig. 3 (a) 
for a standard harmonic D-wave (A((f>) = Ao cos(20)) for 
comparison. The parameters (3 = 1.74 and A /T sc = 2.5 
are used for the harmonic D-wave case. Compared to 
the AFM+D-wave case, we notice the main difference 
in the DOS at low energies shown in the inset. The 
resonance feature of unitary scattering in the residual 
DOS N(u> — ► 0) is seen more clearly but the value of 
N(u> = 0) is smaller than for the AFM+D-wave case 
(see Fig. 3 (a)) with the same impurity scattering rate. 
As in the AFM+D-wave case, l/Ti produces the T- 
linear behavior at low temperatures with impurities and 
the T 3 dependence at higher temperatures below T sc . 
However, there are also some noticeable differences com- 
pared to the AFM+D-wave case. First, the T-linear re- 
gion and the T 3 region arc more clearly separated (a 
smaller crossover region compared to the AFM+D-wave 
gap case). This difference reflects the different shapes of 
the DOS for each case. Second, another more important 
difference is that the magnitude of the low temperature 
l/Ti is about four times smaller for the harmonic D- 
wave than the one of the AFM+D-wave case with the 
same impurity concentration. Therefore, the comparison 
with the data of CeRhIn5 at 1.6 GPa 7 (green diamonds) 
shows a large deviation at low temperatures although 
the T-linear power law is reproduced. It is possible to 
fit the 1.6 GPa data with a much larger impurity scat- 
tering rate of T/Aq ~ 0.16 (not shown in Fig. 3 (b) but 
which can be extrapolated from the shown results). On 
the other hand, in order to fit the 2.1 GPa data (magenta 
diamonds), T/Aq needs to be 0.016 or smaller (see Fig. 3 
(b)). Thus, we would need an explanation for why the im- 
purity scattering rate is suddenly reduced by more than 
one order of magnitude when changing pressure from 1.6 
GPa to 2.1 GPa. 

Summing up the results of our calculations for unitary 
impurity scattering, we found that both gap functions 
with unitary impurities have difficulties to explain con- 
sistently the experiments at 1.6 GPa and 2.1 GPa. 

Now we consider the effects of Born impurity scatter- 
ing. Fig.4 (a) shows the DOS and l/Ti of the AFM+D- 
wave gap with Born limit impurities (c = 1). The inset 



of Fig.4 (a) shows the DOS of the AFM+D-wave gap 
(?7 = 0.3) with different impurity concentrations. Com- 
pared to the unitary impurity case (Fig. 3. (a)), the shape 
of the DOS looks almost the same except for the fact 
that the amount of impurity concentration needs to be 
about five times larger to achieve a similar residual DOS. 
As a result, l/Ti looks similar as in the unitary impurity 
case (Fig. 3 (a)) with five times larger impurity concen- 
trations. We can fit the 1.6 GPa data of CeRhlns with 
the Born impurity scattering rate T/Aq = 0.32. 

However, the results for the harmonic D-wave gap 
with Born impurities are quite different. Fig.4 (b) shows 
the DOS and l/Ti for the harmonic D-wave case with 
Born limit scatterers. With increasing impurity scatter- 
ing rate up to T/Aq = 0.32, the zero frequency residual 
DOS N{0) is less than 20% of N . As a result 1/Tj 
never displays any noticeable T-linear behavior up to 
r/Ao = 0.32. On the other hand, the comparison with 
the 2.1 GPa data of CeRhlns 7,8 shows that the experi- 
mental data can be fitted nicely with T/Aq = 0.16. 

A consistent explanation of the experiments of 
CcRhlns with our calculations is the following. We first 
estimate Ao,2.iGPa/Ao,i.6GPa ~ 2 from the experimental 
observation that T sc -.2.iGPa ~ 2.3K and T S c-,i.6GPa ~ 1 
K. Then at 1.6 GPa where the AFM and SC coexist, the 
SC gap function is the AFM+D-wave gap and we can fit 
the data with Born limit impurities of T/Ao = 0.32 (Fig.4 
(a)). At 2.1 GPa where the AFM disappears, the SC gap 
becomes the standard harmonic D-wave gap. Assuming 
that the Born limit impurity concentration remains the 
same as in the 1.6 GPa sample, the effective impurity 
scattering rate reduces to T/Aq = 0.16, hence the data 
fits well again with our calculation (Fig.4 (b)). 

The remaining question is the origin of the Born limit 
scatterers and its amount. We suggest it could be due to 
strain in the lattice induced with pressure and/or AFM 
domain walls in CeRhlns and chemical substitution in 
CeCu2(Sio.98Geo.o2)2- The scattering rate r/A = 0.32 
in CeRhlns at 1.6 GPa appears quite large. But consid- 
ering Ao oc T sc ~ 1 K, this impurity scattering rate is 
actually quite possible even in the nominally clean sam- 
ple. In fact, the measured specific heat coefficients at low 
temperatures (C(T)/T) both at 1.65 GPa and 2.1 GPa 

ormal 

-0.4 and A 2 .i(0)/A n 

orrnal ~ 0) 14 are 

consistent with our calculations of the residual DOS (see 
the insets of Fig.4 (a) and Fig.4 (b)). 



IV. CONCLUSION 

In summary, we considered the problem of the coex- 
istence of AFM and D-wave SC. Assuming the case of 
A afm > Asc, the gap equation is simplified and the 
AFM correlation imposes an additional constraint on the 
D-wave gap function. As a result we found an interest- 
ing modification of the D-wave gap function; the D-wave 
OP should develop additional nodes besides the original 



4 



D-wave nodes. The DOS of this AFM+D-wave gap solu- 
tion has a generic feature of being almost gapless down 
to very low energy. This shape of the DOS is intrinsically 
vulnerable to a small amount of impurities to produce a 
final DOS having a large amount of gapless excitations. 

We then calculated the nuclear spin-lattice relaxation 
rate 1 jT\ with this AFM+D-wave gap solution as well as 
with a standard harmonic D-wavc gap, and compared the 
results with the experimental \jT x of CeRhIn 5 7 < 8 . We 
found that with unitary impurities both the AFM+D- 
wave gap and a standard harmonic D-wave gap can pro- 
duce T-linear behavior in 1/T\ at low temperature with a 
small amount of unitary impurities. However, both mod- 
els with unitary impurities cannot explain the 1.6 GPa 
and 2.1 GPa data together in a consistent way. Then 
with Born limit impurities, we could fit the CcRhIn.5 1 jT\ 
data at 1.6 GPa with the AFM+D-wave gap with the 
amount of impurities T/Ao = 0.32. Assuming that the 
gap function is changing from the AFM+D-wave form at 
1.6 GPa to the standard D-wave form at 2.1 GPa, we 
could also fit the 1/Ti data at 2.1 GPa with the same 
amount of impurities. The success of our explanation of 
the l/Ti nuclear spin-lattice relaxation rate of CeRhIn 5 
with pressure supports the idea that a general D-wave 
gap function should develop additional nodal points be- 
sides the generic D-wave nodes in the region of coexisting 
AFM and D-wave SC. These additional nodes should be 
observable in carefully designed experiments probing the 
nodal quasiparticlcs 15 . 
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FIG. 1. The schematic 

magnetic Brillouin zone (BZ) (dashed line) and the Fermi 
surface (FS) in the two-dimensional BZ. The crossing points 
of the magnetic BZ and FS are the AFM nodal points. Also 
shown is the schematic gap solution of the AFM+D-wave gap 
equation (The red and blue colors indicate the sign-changing 
OP). 
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FIG. 2. The normalized DOS N(lu)/N for different posi- 
tions of the AFM nodal points 4>afm\ the red dotted line is 
for rj = \4>afm — 4>d\/4>d=0-2) and the blue solid line is for 
rj = 0.4. Inset: The corresponding gap functions A(<^>). 
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FIG. 3. (A) The normalized 1/Ti for AFM+D-wave gap 
(with the AFM nodal point parameter 77=0.3) for uni- 
tary scatterer (c = 0) with impurity scattering rates of 
T/Ao = 0.064,0.032,0.016,0.0. The green diamonds are the 
normalized 1.6 GPa experimental data and the magenta di- 
amonds are for 2.1 GPa data of CeRhIn 5 7 ' 8 . Solid lines 
for T 3 and T are guides for the eyes. Inset: The cor- 
responding normalized DOS N(u))/No, in decreasing order 
of N(ui = 0), with r/A = 0.064,0.032,0.016,0.0. (B) 
The normalized 1/Ti for a standard harmonic D-wave gap 
for unitary scatterer (c = 0) with impurity scattering rates 
T/Ao = 0.064,0.032,0.016,0.0. The green and magenta dia- 
monds are the same data as in Fig. 3 (A). Inset: The corre- 
sponding normalized DOS N(uj)/No. 



FIG. 4. (A) The normalized 1/Ti for AFM+D-wave 
gap (with the AFM nodal point parameter r/=0.3) for 
Born scatterer (c = 1) with impurity scattering rates of 
r/Ao = 0.32,0.16,0.08,0.0. The green diamonds are the 
normalized 1.6 GPa experimental data and the magenta di- 
amonds are for 2.1 GPa data of CeRhIn 5 7 ' 8 . Solid lines 
for T 3 and T are guides for the eyes. Inset: The cor- 
responding normalized DOS N(u))/No, in decreasing order 
of N(lu = 0), with r/A = 0.32,0.16,0.08,0.0. (B) The 
normalized 1/Ti for a standard harmonic D-wave gap with 
Born scatterer (c = 1) with impurity scattering rates of 
r/A = 0.32,0.16,0.08,0.0. The green and magenta dia- 
monds are the same data as in Fig. 4 (A). Inset: The cor- 
responding normalized DOS N(u))/Nq. 
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